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We study the effects of different "fat link" actions for Kogut-Susskind quarks on 
flavor symmetry breaking. Our method is mostly empirical - we compute the pion 
spectrum with different valence quark actions on common sets of sample lattices. 
Different actions are compared, as best we can, at equivalent physical points. We 
find significant reductions in flavor symmetry breaking relative to the conventional 
or to the "link plus staple" actions, with a reasonable cost in computer time. We 
also develop and test a scheme for approximate unitarization of the fat links. While 
our tests have concentrated on the valence quark action, our results will be useful in 
designing simulations with dynamical quarks. 



1 Introduction 



As lattice QCD computations evolve towards quantitative predictions, it becomes increas- 
ingly important to make better use of computational resources by using more sophisticated 
discretizations of the continuum action, or "improved actions" . With Kogut-Susskind quarks 
the errors from the lattice discretization of the fermion action are proportional to a 2 , where 
a is the lattice spacing. In contrast, for the (unimproved) Wilson quark formulation, the 
errors are proportional to a 1 ^. However, for the lattice spacings that are now practical for 
calculations with dynamical quarks, the effects of these errors in the Kogut-Susskind formu- 
lation are not small. In particular, the breaking of flavor symmetry is a large effect. The 
Kogut-Susskind formulation naturally describes four flavors of quarks, which is convention- 
ally reduced to two in dynamical simulations by taking the square root of the determinant. 
However, only a U(l) subgroup of the original SU(4) chiral symmetry is an exact symmetry, 
so for nonzero lattice spacing only one of the sixteen pseudoscalar mesons has a vanishing 
mass when the quark mass vanishes. In principle, flavor symmetry breaking is also present for 
vector mesons, nucleons, etc., but these effects are generally small, and we can concentrate 
on the pseudoscalar mesons. 

Flavor symmetry breaking can be reduced by "fattening" the links, which means that the 
conventional parallel transport using the link matrix is replaced by an average over paths 
connecting the points. This was demonstrated in Ref . Jj] , where the single link was replaced 
by an average of the simple link and three link paths, or "staples". This modification was 
introduced as the simplest possible gauge invariant modification of the nearest neighbor 
coupling, and the relative weight of the staples and simple connection was treated as a 
free parameter. The improvement of flavor symmetry was found to be fairly insensitive to 
the value of the weighting coefficient. Lepage pointed out that this improvement could be 
understood as an introduction of a form factor suppressing the coupling to high momentum 
gluons which scatter quarks from one corner of the Brillouin zone to another [Q]. Lagae 
and Sinclair used this understanding to construct an action where successive smearings 
in all four directions canceled the tree level couplings to all gluons with any momentum 
component equal to n/a, and showed that this further reduced flavor symmetry breaking[||]. 
This action includes paths to nearest neighbor points that are up to nine links in length. 
A fat link Kogut-Susskind fermion action, motivated by a perfect Kogut-Susskind fermion 
action, was tested by Bietenholz and DilgerQ] on the 2D Schwinger Model. They found 
that their action shows very good scaling for the pion and the eta masses. Furthermore, 
perfect action motivated fattening has been used by DeGrand|| for Wilson-like fermions, 
showing very good scaling and renormalization factors very close to unity. In an earlier 
paper we investigated the single staple fattening, together with the Naik term|J improving 
the dispersion relation in more detail [fj]]. Ref. [|7| also describes an algorithm for using these 
actions in full QCD simulations, and extends studies of flavor symmetry breaking to the 
nonlocal pions, which turn out to have much larger breakings than the local non-Goldstone 
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pion which was used in the earlier studies. 



This paper extends our studies of flavor symmetry breaking to actions involving paths 
longer than three links. Many of these results were briefly presented at the Lattice-98 
conference 0. We investigate actions that involve five and seven link paths to the nearest 
neighbor point, in addition to the one and three link paths. The coefficients of the various 
paths are chosen to minimize or to completely eliminate the tree level couplings to gluons 
with transverse momentum components 7t/a. Motivated by encouraging results using links 
fattened by "Ape smearing" |J , we also try a variant of the fattened action which makes the 
fattened links approximately unitary. 



Recently Lepage has analyzed flavor symmetry breaking using the language of Symanzik 
improvement |10|| . In addition to providing a simple construction of the seven-link action 
which cancels tree level couplings at momentum tt, Lepage points out in this work that the 
fattening which reduces flavor symmetry breaking introduces an error of order a 2 in the low 
momentum physics (effectively an a 2 p 2 error, where p is a momentum), and shows how this 
error can be removed by an additional term in the action. 



This work concentrates on the effect of the valence quark action on the flavor symme- 
try breaking. We have not done a study of the effects of changing the dynamical quark 
action; however, we have found in previous studies that improvements in flavor symmetry 
arising from changes in the valence quark action are indicative of improvements in full QCD 
simulations^. Also, we have not studied scaling, or the independence of mass ratios on the 
lattice spacing, although we expect that they will show improved scaling properties. Our 
concentration on the splittings of the pions is motivated by experience that suggests that 
this is the worst practical problem with Kogut-Susskind quarks. (For example, in Ref. |7j we 
found that rotational symmetry of the pion propagator was essentially restored at coarser 
lattice spacings than was the flavor symmetry.) Furthermore, we expect that the study of 
the pion spectrum provides a good guide to the quality of actions, since the low energy 
dynamics of QCD is approximately the dynamics of a pion gas. 



2 The actions 



Figure [l] illustrates the coupling to gluons with momentum components ±7r. The figure 
illustrates the single link couplings in D^, connecting the central point to both the forwards 
and backwards nearest neighbors, as well as three link staples connecting to the same points. 
Each horizontal arrow represents a parallel transport by ~ 1 + igaA^. The directions of 
the arrows in the coupling to the backwards neighbor include one minus sign because this 
transport involves W and another minus sign appearing explicitly in the derivative. Now 
consider the coupling to a //-direction gluon (A^ne-v)) with momentum component ±7r in 
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Figure 1: A simple link and staples in one plane, showing the relative signs of the coupling 
to a gluon with momentum component 7r in the transverse direction (unparenthesized) , and 
a gluon with momentum 7r in the longitudinal direction (parenthesized). 



the v direction. The ± signs above the /x direction links indicate the sign of the coupling 
of this gluon relative to the forward simple link. Thus, if C\ and C3 are the weights of the 
one-link and three-link paths, the coupling to this gluon from the paths pictured here is 
Ci — 2C3. (This is not the whole story — one must also include the staples in the directions 
that are not shown in this figure.) Similarly, the parenthesized signs below the /x direction 
links show the relative coupling to a gluon with momentum 7r in the /x direction (A^ire^)). 
Note that this coupling is automatically cancelled between the forward and backward parts 
of Dp, so we don't have to worry about longitudinal momentum 7r. Now one might worry 
that the v direction links, required to keep the expression gauge invariant, might introduce 
couplings to v direction gluons with momentum components ±7r. But for A v {i:e^) gluons, 
the contributions of the vertical links in the center cancel, and the contribution of the links 
from the left and right sides cancel, since they are separated by 2a and traversed in opposite 
directions. Similarly, the coupling to A u (ire u ) gluons cancels between the top and bottom 
halves of the figure. This argument extends to the five and seven link paths illustrated in 
Fig. ^, so the end result is that we can compute the couplings to gluons with momentum 
components equal to 7r by just considering the /x direction links in D^. Taking into account 
the multiplicities of the various paths and using C n for the weight of the n link paths, the 
couplings are given below. In these expressions the coefficients explicitly show how many 
paths of each length contribute with positive weight and how many with negative weight. 
We also write explicitly (2C5) and (6C7) to indicate that in Jp x there are two shortest paths 
connecting the starting point to the x direction link displaced by +y + z, and six paths to 
the link displaced by +y + z + 1. 



• Coupling to k = (0, 0, 0, 0): d + 6C 3 + 12(2C 5 ) + 8(6C 7 ) 
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Figure 2: The simple link, three link staple, five link staple and seven link staple used in 
suppressing flavor symmetry breaking. The final path is the five link path used, following 
Lepage, to correct the small momentum form factor. 

• Coupling to k = (0, 7T, 0, 0): C x + (4 - 2)C 3 + (4 - 8)(2C S ) + (— 8)(6CV) 

• Coupling to k = (0, 7T, 7T, 0): Cx + (2 - 4)C 3 + (4 - 8)(2C 5 ) + (+8)(6C 7 ) 

• Coupling to As = (0, vr, vr, tt): d + (-6)C 3 + (12)(2C 5 ) + (-8)(6C 7 ) 

If we are willing to use all the paths up to length seven, we can normalize the zero momen- 
tum coupling to one and set all the others to zero with G\ = 1/8, C 3 = 1/16, (2C 5 ) = 1/32 
and (6C7) = 1/64. This defines our "Fat 7" action and, with tadpole improved coefficients, 
the "Fat7tad" action. However, it is interesting to ask how helpful this complexity is - 
could we get by with shorter paths? If we restrict ourselves to five link and shorter paths, we 
can no longer satisfy all of these conditions. However, we can choose the couplings to mini- 
mize the maximum of the couplings to the high momentum gluons. This leads to Ci = 1/7, 
C 3 = 1/14 and (2C 5 ) = 1/28, with all couplings to high momentum gluons reduced by a 
factor of seven. This defines our "Fat5" action. 

Here we note that in two dimensions the equivalent of our "Fat 7" action has only a 3-link 
staple. Our tree level formula for the relative weight of the staple would be 0.25 in two 
dimensions. This is very close to the relative weight of 0.238 introduced in the approximate 
perfect action constructed in Ref. [Q for the 2D Schwinger Model. 

All of these actions can be tadpole improved by inserting a factor of (1/mo) L1 i n the 
coefficient of each path, where L is the length of the path. Using L — 1 instead of L amounts 
to absorbing one power of u into the quark mass. We use the average plaquette to define 

Uq. 
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While the paths shown in Fig. ^ can be used to reduce or eliminate couplings to gluons 
with transverse momentum components equal to 7r, Lepage has pointed out that they have 
the undesirable effect of modifying the coupling to gluons with small nonzero transverse 
momentum components, essentially by introducing a second derivative coupling proportional 
to a 2 p 2 |l(]]. Following Ref. [[Kj, we can correct for this by introducing a flavor conserving 
five link path, +y + y + x — y — y, into D x , giving an action with no tree level order a 2 
corrections. While it is not clear what observable quantities are affected by this correction, 
it is a relatively cheap and aesthetically pleasing addition to the Fat7+Naik action. We have 
tested flavor symmetry breaking with this action, which we call the "Asq" action, or, when 
the coefficients are tadpole improved, the "Asqtad" action. 

To summarize and clarify these actions, we give the coefficients of the paths in the a 2 
improved action in a form useful for simulation. Here c\ is the coefficient of the simple link. 
C3, C5 and C7 are the coefficients of the three, five and seven link paths in Fig. Q is the 
coefficient of the three link path to the third nearest neighbor (Naik term), and cl is the 
coefficient of the five link path implementing the correction introduced by Lepage. The origin 
of each term in the coefficients is identified by a subscrpt: "F" for flavor symmetry, "N" for 
the Naik correction to the quark dispersion relation, and "L" for the small momentum form 
factor correction. The factors of 1/2 and 1/6 in C5 and Cj compensate for the number of 
different paths connecting the starting point to the /z direction link parallel to the simple 
link. For example, the x direction link displaced by +y + z (coefficient C5) is included in 
both +y + z + x — z — y and +z + y + x, — y — z, and we average over both paths. 

Cl = (l/8) F +(3/8) L + (l/8) w 

c 3 = (1/16) F 

c 5 = (1/32)^(1/2) 

c 7 = (l/64) F (l/6) 

c L = (-1/16) L 

c N = (-1/24)* (1) 



We have also experimented with an action in which the fat links are approximately 
unitary. Theoretically, it is not clear why unitarity should be a concern in suppressing flavor 
symmetry breaking. However, results in Ref. || led us to consider "APE-smeared" links, 
where an average over the single link and nearby paths is constructed as above, and then the 
resulting fat link is projected back onto the nearest element of SU(3). Of course, such actions 
are not explicitly expressed as a sum over paths, and dynamical simulations using them would 
require an extension of our algorithms. However, spectrum calculations using APE-smeared 
links for the valence quarks are straightforward. In particular, we ran a calculation using 
a single iteration of APE smearing on a fairly coarse lattice. This action differs from the 
one link plus staple action only in the projection onto a special unitary matrix, yet, as will 
be seen in a later section, it produced smaller flavor symmetry breaking than the one link 
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plus staple plus Naik action. (The Naik term has little effect on flavor symmetry.) This 
motivated us to construct an action for which the fat links are approximately unitary, but 
are still expressed as an explicit sum over paths. 



An almost unitary matrix M can be expressed as a unitary matrix U times a correction: 

M = U(l + e) , (2) 
where e is Hermitian. Now look at M^M to first order in e 

M f M = (1 + e)U^U(l + e) w 1 + 2e . (3) 
Now, inverting Eq. |2| to first order, 

U = M (1 -e) = -M - l -MM^M . (4) 

We want to use this equation to make a fat link approximately unitary. Let M be a generic 
fat link: 

M = L + aS (5) 

where L is the simple link and S is some sum over other paths connecting the ends of L. 
S may be the sum of the staples for a simple fattening, or something more complicated 
involving longer paths. We will work to first order in a. Customarily we also rescale M by 
some factor which is 1 + ba. However, this factor cancels in the derivation, so we suppress 
it here. Unitarize M as above to make an approximately unitary fat link F: 

F = -M — -MM l M 

2 2 

3 (L + aS) - I (L + aS) (tf + aS^) (L + aS) (6) 



2 y '2 

Expand, keeping terms up to order a and using L^L = 1, and remarkable cancellations occur 

F = L+Z(S-L&L) (7) 



F is an approximately unitary fat link expressed as a sum over paths — a form suitable 
for dynamical simulations. The only change is that the sum over fattening paths, S, has 
been replaced by an average of the paths traversed in each direction, | (S — , with L 
inserted as necessary to maintain gauge invariance. At this point one can forget about how 
it was derived, and just verify from Eq. |7| that F^F = 1 to order a. It is also easy to verify, 
by setting link matrices to 1 + iA^, that F averages over position in the same way that 
M did. The minus sign in front of the LS^L term compensates for the fact that this path 
is now traversed in the opposite direction. This action is illustrated in Fig. ||| for the case 
where S is the three link staple. 
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Figure 3: An approximately unitary variant of the three link staple. 

The two terms can be made to look more symmetric by factoring L out on the right (or 
left). 

F = (l + £ (Stf - L£t)) L (8) 

In this form, we are just following the parallel transport by L with a difference of going 
around closed loops in opposite directions, or a term proportional to F^ u . 



3 The simulations 



Most of our spectrum calculations used two of the same sets of sample lattices used in Ref. J7|. 
These were lattices with two flavors of dynamical quarks, where the dynamical quarks used 
the " Staple+Naik" action, and a Symanzik improved gauge action at 10/gf mp = 7.3 and 7.5. 
At 10 /g? = 7.3 we used 12 3 x32 lattices with dynamical quark masses of 0.02 and 0.04, with 



lattice spacings determined from tq [11 of about 0.15 fm. and 0.16 fm. respectively. The 



10/gf mp = 7.5 runs used 16 x 48 lattices with m q = 0.015 and 0.030, with lattice spacings 
of 0.13 fm and 0.14 fm, respectively. Sample sizes ranged from 48 to 60 lattices]?]], with four 
source time slices per lattice for spectrum calculations. In addition we computed the meson 
spectrum on a set of eleven large (32 3 x 64) quenched lattices with a lattice spacing around 
0.07 fm., using the one plaquette gauge action at Q/glonv = 6.15. 

For each of the three gauge couplings we computed spectra at two values of the quark 
mass, which allows us to interpolate results to make fair comparisons of the actions. (For the 
two flavor runs, the sea quark mass was also changed.) The meson masses were interpolated 
assuming that the squared meson masses are linear in the quark mass. One approach to 
comparing actions is to interpolate to a fixed ratio of nic to m p , where mc is the Goldstone 
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pion mass. In so doing, we are letting each action "determine its own lattice spacing" . Since 
the p mass is also dependent on the valence quark action, this results in a slightly different 
estimate of the lattice spacing for each action. An alternate approach is to assume that each 
set of lattices has a fixed lattice spacing, which could be determined either from the rho mass 
for a fixed choice of action, or from the static quark potential. Operationally, this results 
in interpolating all the valence actions on a set of lattices to the same pion mass, and we 
choose the Goldstone pion for this interpolation. For the runs with lO/gf mp = 7.3 and 7.5 we 
will present results with both sets of assumptions. For the fine lattices, with 6/gl onv = 6.15, 
the differences in the p masses with different actions are insignificant, which is expected as 
the lattice spacing gets smaller. 

As in Ref. 0, we parameterize the splitting of the pions by the dimensionless quantity 

2 2 

ml — Vfln , , 

h = — 2 f , 9 

m p ~ m G 

where m w is one of the non-Goldstone pion masses, mc is the Goldstone pion mass, and 
m p is one of the local p masses. Since the p masses are nearly degenerate, it makes little 
difference which one we use. In our analysis we used the local ji (g> ji (p 2 ) mass because for 
heavy quarks it is often estimated more accurately. We use the squared meson masses since 
they are approximately linear in the quark mass and we interpolate and extrapolate in quark 
mass. Empirically 8% is fairly insensitive to the quark mass, and the theoretical analysis in 
Ref. E2| predicts that the numerator is independent of quark mass for small quark masses. 



The squared Goldstone pion mass is subtracted in the denominator to give more sensible 
behavior at large quark mass. At large am q , all of the meson masses become degenerate, but 
it is still possible to use the difference between the Goldstone pion mass and the rho mass 
as a natural scale for flavor symmetry breaking. 



4 Results 



Figure |] shows the spectrum of pion masses obtained with the "Staple+Naik" action at 
10/ gf = 7.3 and 7.5 respectively. Here, as in all of our calculations, the pattern of near- 
degeneracies predicted in Ref. |l2j is evident. Note that the local non-Goldstone pion, the 
flavor 7q75 pion, is one of the lightest non-Goldstone pions, so a realistic assesment of the 
flavor symmetry breaking requires consideration of the nonlocal pions. All to the actions we 
have tested produce a qualitatively similar pattern of pion masses. 
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Table 1: Path coefficients for the various quark actions. The first column is a name for 
the action. The remaining columns are the coefficients for the paths in Fig. Q and for the 
connection to the site at distance three (Naik term). Notice that in the tadpole improved 
actions we have absorbed one power of uq into the quark mass, so uq appears to the power 
length — 1. 

Lepage Naik 



9 ^1 
8 27 



Action 


Link 


3-staple 


5-staple 


7-stapl 


OneLink 


1 








Staple+Naik 


9 1 

8 4 


9 1 

8 8 






Fat5 


1 
7 


1 1 

7 2 


i i 

7 8 




Fat5tad 


1 

7 


I I,,- 2 


7 8 M 




Fat7 


1 

8 


1 1 

8 2 


1 1 

8 8 


1 1 

8 48 


Fat7tad 


1 

8 


8 2^0 


8 8^0 


8 48 M 


Asqtad 


8 ~ 8 ~ 8 


8 2 U 


I I v/ - 4 

8 8 U 


8 48 "0 



16 "0 24 "0 



Figure || shows 52 from three different actions, each with three different lattice spacings. 
All of these results are interpolated in quark mass to the point where m^/mp = 0.55. The 
improvement of the flavor symmetry breaking is evidenced by the smaller non-Goldstone 
pion masses for the improved actions. In this figure we have determined the lattice spacing 
for each action using the p mass evaluated with this action. Thus, even though the spectra 
for the different actions are evaluated on the same sets of lattices, they appear at different 
horizontal positions. However, as we would hope, the ambiguity in lattice spacing coming 
from the difference in p mass with action becomes smaller at the smaller lattice spacings. (It 
is neither surprising nor upsetting that the p mass should depend on the action. After all, 
in the end we expect that these actions will have better scaling behavior in all quantities, 
meaning that we should come closer to the continuum limit at these large lattice spacings. 
Among other things, this means that we expect the nucleon to rho mass ratio, a well known 
problem on coarse lattices, to be improved with these actions.) 

In Table |] we show 82 for the various actions. We tabulate this for only three of the 
pions: the local non-Goldstone pion with flavor structure 7075, and the two three-link pions, 
with flavor structures 70 and 1. The local pion is the most often studied one, and so allows 
comparison with other work. The flavor 1 is the worst case pion. The other three link pion 
is in the second worst multiplet, but it is interesting to study because its parity partner 
would have exotic quantum numbers, so the propagator can be fit to a simple exponential, 
leading to smaller errors for the mass estimates. We note in passing that the flavor 1 pion 
is properly called a "pion" here instead of an "eta", since we did not compute quark-line 
disconnected diagrams in the propagator. In other words, one may imagine that the quark 
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and antiquark in this pion carry a flavor quantum number in addition to that coming from 
the Kogut-Susskind quarks' natural four flavors, and so cannot annihilate each other. Thus, 
it is correct to demand that an improved action should make this pion degenerate with 
the others. In order to make a fair comparison of the actions, we have interpolated (or 
extrapolated) the spectrum to a fixed mc/wip = 0.55 point. This is done using two different 
bare masses for the valence quarks. We also interpolated (or extrapolated) our spectrum 
to a fixed me = 545MeV using the heavy quark potential parameter r in order to fix the 



scale |jTlj]. As we can see from the table, the computed <Vs in both cases are in agreement 
with in errors. The errors in both cases were computed with jackknife analysis. Since 8 2 
involves a ratio of mass differences, and the different masses are all correlated, naive error 
propagation would lead to an overestimate of the errors in 8 2 - Therefore we used a jackknife 
analysis to compute these errors. When comparing actions, we are interested in the difference 
in 82 between the two actions. Since we used the same lattices for all the actions, these 82 's 
are not independent, and one should really do a jackknife analysis in order to determine the 
errors on the differences of 8 2 , s. We have done this and we find in general the same error as 
the one computed with error propagation from the quoted 82 errors. 



Our results show that all the variants of the fat actions significantly improve flavor 
symmetry, with larger improvement for the actions that suppress more couplings to gluons 
with transverse momentum tt. Generally, tadpole improving the tree level coefficients results 
in a better action. The approximately reunitarized action with coefficient 0.25 seems to 
work better than the Staple+Naik action. We have no clear understanding of why this 
happens, but this may be a hint that a more careful (non-perturbative) tuning of the fattening 
coefficients would result in better actions. As expected, the Fat7tad action is the best for 
suppressing flavor symmetry violations. The Asqtad action is slightly worse than the Fat7tad 
action in this regard. Extra 0(a 4 ) flavor symmetry violation introduced by the Naik term 
and the Lepage term may be responsible for this. However, the Asqtad action has improved 
rotational symmetry due to the Naik term, and has no additional 0(a 2 ) flavor conserving 
errors introduced by the fattening. Thus, of all the actions we have studied, we consider this 
one to be the best candidate for an improved Kogut-Susskind action. 

In Ref. m it was shown that smearing the links by 'Ape smearing", where the link 
and staples are averaged and the result in projected back onto SU(3) (i.e. replaced by the 
SU(3) matrix which maximizes Tr(WF), where F is the fattened link) improved the flavor 
symmetry, among other nice features. We find similar results on our set of sample lattices. 
It is interesting to compare the "Apel" action in Table ^ with the "Staple-un(.25)" action, 
since the difference between the two is that the Apel action uses links that are exactly 
unitary, while the Staple-un(.25) uses the same fattening, but is only unitary to first order 
in the staple weight. We see that the approximately unitary action is only slightly worse 
than the Apel action. Also, we see as expected that that the Ape4 valence action gives 
a very good suppression of flavor symmetry breaking, although using it for a dynamical 
action is difficult. The approximately reunitarized actions seem to have slightly better flavor 
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symmetry breaking than the comparable Link+Staple action, but probably not enough to 
justify the extra complexity 

Another practical advantage of the fat link actions is that the fat link configurations 
are smoother than the original configurations, and the conjugate gradient computation of 
the propagators converges in fewer iterations. Quantifying this statement is tricky because, 
unlike the flavor symmetry breaking, where 82 is approximately independent of quark mass, 
the number of conjugate gradient iterations is very sensitive to the quark mass. For ex- 
ample, at lO/gf mp at a fixed bare quark mass of 0.02 the Link+Staple action requires 21% 
fewer conjugate gradient iterations than the conventional action. However, after interpola- 
tion to nic/rap = 0.55 this advantage disappears, and both of these actions require about 
the same number of iterations. On the coarse lattice the approximately unitary actions do 
better — the Staple-un(.25) action and the Apel action require about 10% fewer iterations 
than the conventional or Link+Staple actions, and the Ape4 21% fewer. On the fine lat- 
tice, with 6/g^ onv = 6.15, there is much less ambiguity since the various pion masses are 
much closer. Also, the fattening does a much better job of smoothing the configurations on 
the fine lattice. In this case the Link+Staple action requires 46% fewer iterations than the 
conventional action, and the Fat7 action 52% fewer. However, the advantage of the approx- 
imately reunitarized action has disappeared - it requires just about as many iterations as 
the Link+Staple. Regrettably, the Asqtad action does not do as well as the Fat7 action, 
only reducing the number of iterations by 37% as compared to the conventional action. This 
is probably because of the negative coefficient associated with the Lepage term (see Table 
[3]), meaning that this term is actually undoing part of the smoothing accomplished by the 
other paths. We should note that with these complicated actions the cost of simulation with 
dynamical quarks is no longer completely dominated by the conjugate gradient. Except for 
very light quark masses, the cost of computing the fermion force and precomputing the fat 
links becomes comparable to the conjugate gradient cost. 



5 Conclusions 



In this paper we have investigated the possibility of constructing a Kogut-Susskind action 
with improved flavor and rotational symmetry suitable for dynamical fermion simulations. 
For this reason we want to keep the amount of new paths introduced into the action as small 
as possible. An action containing fat links with paths up to length seven was constructed. 
At tree level this action has no couplings of quarks to gluons with a transverse momentum 
component 7T. As a result, at tree level the flavor symmetry violating terms in the action are 
completely removed. In addition to flavor symmetry, the rotational symmetry is improved 
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Table 2: The flavor symmetry breaking measure 82 for the various actions and lattice 
spacings. Results are shown for the local non-Goldstone (7075) pion and for the two three- 
link pions, with Kogut-Susskind flavor structures 70 and 1. The actions are those specified 
in Table [], plus actions with a single staple fattening approximately unitarized as discussed 
in section 2. For the actions with single staple fattening we give the relative weight of each 
staple to the single link in parentheses (e.g. Staple-un(.25)). The "Apel" action is the single 
link plus the three link staple with a relative weight of 0.25, projected back onto an element 
of SU(3) to numerical precision. For the "Ape4" action, this smearing and projection was 
repeated four times. Otherwise the coefficients of the various paths are found in Table [| 
For comparing the actions we interpolate (or extrapolate) to mc/mp = 0.55 or alternatively 
to m G = 5A5MeV. 
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by introducing the Naik term. Finally, as Lepage pointed out|K|, we need to introduce an 
extra five link staple in order to cancel errors of 0(a 2 p 2 ) introduced by the fattening. The 
resulting action can be further improved by tadpole improvment. This action, which we call 
"Asqtad", is an order 0(a 4 ,a 2 g 2 ) accurate fermion action. 



Asqtad is an action simple enough to be useful for dynamical simulations. Preliminary 
tests using our generic code|7| show about a factor of 4 higher cost than the code implement- 
ing the standard Kogut-Susskind action. We have found optimizations specific to the Asqtad 
action that could bring the cost factor down to 2-2.5. The 0(a 2 ) precision at a cost of a 
factor of 2-2.5 makes such an action very competitive with the other popular improved ac- 
tions such as D234|l3|], perfect or approximately perfect actions, the Neuberger action, and 
domain wall fermion actions. The highly improved chiral symmetry of actions respecting 
(or approximately respecting) the Ginsparg- Wilson relation is not something that Asqtad 
can compete with. On the other hand, cost may favor the Asqtad action. The Neuberger 
action j~5[ , approximately perfect actions [16], |9|, |J and domain wall fermions||17|| are all 



fairly costly to implement for dynamical fermions. 



In view of the enormous price one has to pay in order to have highly improved chiral 
symmetry on the lattice, we think that the Asqtad action is a good candidate for a fermion 
action to be used in the next generation of dynamical simulations. Flavor symmetry breaking 
in the traditional Kogut-Suskind action at lattice spacings commonly used in high tempera- 
ture QCD studies results in pions as heavy as the kaons, making it impossible to study the 
effects of the strange quark. Our study of the Asqtad action shows that one can achieve a 
good separation between the pions and the kaons at accessible lattice spacings. Thus Asqtad 
is an action that may prove very useful in projects in which the effects of the strange quark 
are to be studied. 
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Figure 4: The spectrum of pion masses with the "Link+Staple+Naik" action. The two 
spectra are at 10 /gf = 7.3 and 7.5 respectively. In both cases we have interpolated in 
quark mass, assuming that squared meson masses are linear in the quark mass, to the point 
where mc/ra p = 0.55. The pion masses are plotted in units of the p mass. Each pion is 
labeled by the gamma matrix specifying its flavor structure. The pattern of near degeneracies 
predicted in Ref. [f2] is evident. As expected, the mass splittings among the pions are smaller 



for the smaller lattice spacing (10/ gf = 7.5) 
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Figure 5: The 5 2 for three different actions versus the squared lattice spacing. From left to 
right, the actions are the conventional "One Link" action, the "Link+Staple+Naik" action, 
and the "Order a 2 tadpole" action. In each case, the p mass for the respective action was 
used to define the length scale. The points come from two flavor runs with 10 /gf = 7.3 
(rightmost points in each graph), 10/gf mp = 7.5, and a quenched simulation with single 
plaquette gauge action at Q/g1 onv = 6.15. The lines drawn in each panel are approximate 
slopes for the local non-Goldstone pion (lower line) and for the three link flavor 1 pion (upper 
line), which is the worst case pion. 
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